Elasticity model of a supercoiled DNA molecule 
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Within a simple elastic theory, we study the elongation versus force characteristics of a supercoiled 
DNA molecule at thermal equilibrium in the regime of small supercoiling. The partition function is 
mapped to the path integral representation for a quantum charged particle in the field of a magnetic 
monopole with unquantized charge. We show that the theory is singular in the continuum limit 
and must be regularised at an intermediate length scale. We find good agreement with existing 
experimental data, and point out how to measure the twist rigidity accurately. LPTENS 97/28. 



The measurements on single DNA molecules, beside 
their possible biological interest, provide a wonderful lab- 
oratory for the physical studies of a single polymer chain. 
For instance, recent experiments have shown that the 
elongation versus force characteristics of a single DNA 
molecule [^| is very well fitted ||] by the well known 
worm-like chain (WLC) j| which describes a chain by 
an elastic continuous curve at thermal equilibrium, with 
a single elastic constant, the persistence length A charac- 
terizing the bending energy. The WLC can be solved an- 
alytically by mapping it to a quantum mechanical prob- 
lem. Its partition function is nothing but a Euclidean 
path integral for a quantum dumbbell, which can be com- 
puted, in the relevant limit of long chains, by finding the 
ground state of the corresponding Hamiltonian. 

Our work is motivated by the more recent experiments 
which have measured the elongation versus force charac- 
teristics of a supercoiled DNA molecule || . We use the 
simplest generalization of the WLC with twist rigidity, 
which is supposed to work at small supercoiling angles. 
The DNA molecule is described by a thin elastic rod in- 
volving a new elastic constant, the twist rigidity C. The 
stress imposed by twisting the endpoint of an open rod 
can be absorbed both in some twist (if the axis of the 
rod is a straight line), and in some deformation of the 
rod's axis. This decomposition is well known in the case 
of closed rods, for which the experimental constraint is 
a topological invariant It results in a subtle com- 
petition involving the creation of plectonemes which has 
received quite a lot of attention, both for the study of 
the ground state, and also taking care phenomenolog- 
ically of thermal fluctuations around some low energy 
configurations @,^,|). In contrast, we keep here to the 
simplest regime of small supercoiling, but we provide a 
full analytic and numerical study of the twisted open rod 
at thermal equilibrium at a finite temperature, extending 
thus the standard WLC analysis to this case. While it has 
been known for some time that the elastic thin rod model 
agrees qualitatively with experimental results on DNA, 
this knowledge was based only on some Monte Carlo sim- 
ulations flj,[l3| . Our analytic study of the worm-like rod 
chain (WLRC) shows very good quantitative agreement 
with the experimental data at small supercoiling and al- 



lows for a precise determination of the ratio C/A, which 
was poorly known so far. 

We shall show that beside the two elastic constants 
describing the bending and twisting rigidity, one needs 
to introduce an intermediate lengthscale b between the 
microscopic interbasepair distance and the persistence 
length, which plays the role of a short distance cutoff. 
Obviously the rod description cannot hold below a cut- 
off length scale b which is at least of the order of the 
double helix pitch p ~ 0.34 nm, since our description av- 
erages over these periodic oscillations. As we shall see the 
WLRC continuous limit is singular and shows properties 
qualitatively very different from any discretized version of 
the chain. So the existence of a cutoff, which is expected 
on physical ground, is crucial, but the final results turn 
out to be fairly independent of its precise value, within a 
reasonable range around .5 nm. Similar singularities of 
the continuous limit are well known in the winding prop- 
erties of pure random walks ficj] . Their appearance here 
is not fortuitous since the worm-like rod chain is related 
to random walks in rotation space. It is interesting that 
supercoiled DNA presents an experimental system where 
these subtleties of the continuous limit of random walks 
turn out to be relevant. 

The WLRC, already studied in (l2[|^Q, is described 
in the continuous limit by the orthonormal triedron 
{t(s), u(s), n(s)} where s is the arc length along the 
molecule, t is the unit vector tangent to the chain, and 
n describes the orientation of the rod. For describ- 
ing DNA, this triedron is obtained by applying a rota- 
tion 1Z(s) to a reference triedron which characterizes the 
natural helical structure of the molecule. The rotation 
TZ(s) is parametrized by the usual three Euler angles 
9(s),<fi(s),ip(s), and the reference triedron is such that 
8(s) = 0, 4>(s) + ip(s) = ujqs, where loq is the rotation 
per unit length of the base axis in a relaxed rectilinear 
DNA molecule. With the above definition, the set of 
s dependent Euler angles 8(s),(j)(s),i(j(s) describes the 
general deformations of the DNA molecule with respect 
to the relaxed rectilinear configuration. We neglect self- 
avoidance, the energy of a chain of length L is purely 
elastic: E e i = ds(eb + e*), with the bend and twist 
energy densities given by: 
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where the dot stands for the s derivative. We work in 
units where the temperature ksT = 1, so that A and 
C have dimension of a length. The discretised version 
is defined by quantifying s as an integer multiple of an 
elementary length scale b, and approximating integrals 
and derivatives by sums and differences, while keeping 
the periodicity. We study the equilibrium properties of 
such a rod pulled by a force F = Fx. The total energy is 
thus E — E el — F J L ds cos#(s). 

The partition function of the elastic chain described 
by eq.([T]) is nothing but the Euclidean path integral for 
a quantum symmetric top, with the important difference 
that the eigenfunctions are not periodic in the angles ip 
and 4>. Therefore the momenta conjugate to these angles 
will not be quantized. In our analytical work, we suppose 
for simplicity that the boundary values of the Euler an- 
gles are 0(0) = 6{L) = 0, and we define 0(0) = ip(0) = 0. 
Then the experimentally imposed supercoiling angle \ 
amounts to fixing: ip(L) + <j>{L) = X- In the open DNA 
chain the continuous angle \ replaces the topological 
linking number ■ We shall limit ourselves to configu- 
rations where the Euler angles are regular enough, such 
that x can be written as the integral ds (^ip + <fj . 
It is convenient to introduce the rod twist T w , given by 
the integral T w = ds (ip + <j> cos 9^j ; T w appears as a 
Gaussian variable in the partition function path integral. 
We now define a 'local writhe' contribution to the super- 
coiling angle xw as xw = X — T w = f Q L ds 0(1 — cos 9). 
The above decomposition || is reminiscent of the one 
performed in the case of of closed chain where the topo- 
logical linking number Lk is decomposed into twist and 
writhe @. 

The partition function for a fixed value of x is given 
by the path integral in the space of Euler angles: 

Z = J d[cos 9, 0, V] 5\x - £ ds(j> + j,) J e- E (2) 

After introducing an integral representation of the S func- 
tion which fixes x, one can perform the gaussian path 
integral on the angle i/j. Z is then expressed as a path 
integral on the two angles cos 9 and 0, with an effective 
energy: 
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This form (^) is useful for numerical simulations |Q af- 
ter a proper discretization, but not for analytic computa- 
tion, due to its non local character. Alternatively we can 
compute the x Fourier transform Z = J dx i?exp(— ikx), 



which is again given by a path integral on the two angles 
cos 9 and 0, with the effective energy: 
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This last form has an appealing quantum mechanical in- 
terpretation: If one analytically continues the s-integral 
towards the imaginary axis, one recognizes the action 
integral of a particle with unit charge moving on the 
unit sphere under the joint action of the electric field 
F and the magnetic field A^ = k (1 — cos6>) of a mag- 
netic monopole of charge k. One easily deduces the cor- 
responding Hamiltonian iJ, by substituting — —ijfj; 
by P4, — A<f, in the WLC Hamiltonian (which corresponds 
to Aj, — 0). Because of the averaging over the final 
<p — 4>(L), only the eigenvalue m = of p^ contributes 
and we can set p^ = in H. We work with the dimen- 
sionless quantities H — H/A and a = AF, in terms of 
which we get: 
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Introducing the eigenstates and the eigenvalues of H , 
H^ n (k,9) = e n {k 2 , a)^ n (k, 9), the Fourier transformed 
partition function Z can be written as the sum: 

Z = J2\^n(k 1 0)\ 2 exp(-^ (e n (k 2 ,a) + ^-)) (6) 



In the large L limit, the sum over the eigenstates is 
dominated by the one with lowest energy eo(a,k 2 ), if 
L/A ^> Ae where Ae is the energy gap between the 
ground state and the nearest excited state of H . This 
gives the approximate expression for the partition func- 
tion Z : 
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Therefore one can deduce, from the ground state energy 
eo (a, k 2 ) of the Hamiltonian H, the observable properties 
of a long WLRC, of which we now discuss two important 
ones. The relative extension of the chain in the direction 



of the force is given by (z)/L = {AjL)^§£-. If instead 
of constraining x °ne measures its thermal fluctuations, 
their probability distribution is just P{x) oc Z . For in- 
stance the second moment is given by: 
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This expression shows that the WLRC is pathological 
because of "giant" fluctuations of xw- The contribu- 
tion to < x 2 > from the twist fluctuations, ^, scales 
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linearly in L, as one expects in a one dimension statis- 
tical mechanics system with a finite correlation length. 
In contrast the second piece of (||) giving the contri- 
bution from (xw)> i s divergent: evaluating eo(k 2 ,a) at 
small k 2 from standard perturbation theory we find 
ixw 2 ) = (L/A){ (1 -cos 0)/(l+cos0)) o where ( } is the 
quantum average taken on the groundstate $o {&) of the 
WLC Hamiltonian (which is H at k — 0). As $ (tt) ^ 
(for any finite force), we get a logarithmically divergent 
result. One can show that e(a,k 2 ) ~ e(a, 0) + |fc|$o( 7r ) 2 - 
This linear behaviour of the energy in |fc| shows that P(x) 
has a Cauchy tail, and thus a diverging second moment. 
The Cauchy distribution has been verified in the limit of 
a vanishing force, a = 0, where the eigenfunctions can be 
found exactly in terms of Jacobi polynomials. A related 
consequence is that the extensive part of the average ex- 
tension is unchanged by the supercoiling angle x at small 
forces: < z > /L ~ 2a/3 independently of F, in striking 
contradiction to experiment. 

In contrast to the WLC, the continuous limit of the 
WLRC is singular. This singular behaviour could have 
been anticipated since H describes the motion of a 
charged particle in a magnetic monopole with an un- 
quantized magnetic charge, a notoriously ill defined prob- 
lem if no cutoff is provided near 9 — n. As we have 
argued previously the WLRC model is not expected to 
give a good description of supercoiled DNA unless one 
introduces the cut-off length scale b. The non trivial fact 
is that this existence of a cutoff affects the 'macroscopic' 
properties taking place on the length scale of the whole 
molecule. 

In order to validate the discretized WLRC, we have 
performed a Monte Carlo simulation, mostly using the 
discretized version of (^|). Such simulations are known 
to account well for the observed behaviour of circular 
DNA @, and have been used recently for the study 
of chains elongated with large supercoiling angles [[l3[ . 
With respect to these works, we have discarded the self 
avoidance, since our aim is to test the discretized WLRC 
without self avoidance. We have discretized the chain 
with elementary rods of length b = A/10, and simulated 
mostly chains of length L = 30A In order to facilitate 
the thermalisation, we have relaxed in the simulation the 
constraint 9(L) = 0, which should not affect the exten- 
sive quantities. The elementary moves which we used 
was to scan sequentially each point i = 1, N = L/b in 
the chain, and propose a global rotation of the tangent 
vectors tj,j = i,...,N around a random axis with an an- 
gle 7 taken with a flat distribution in [—70, 70], where 70 
is chosen such that the acceptance rate of the moves is of 
order .5. These rotations of a fraction of the chain were 
the best we found for insuring a relatively fast thermal- 
ization jl4| . The results presented in Fig.l, obtained with 
C/A = 1.4, show that the elongation versus x character- 
istics reproduces well the experimental values at small 
enough x- 
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FIG. 1. The elongation versus reduced supercoiling angle 
7] = X-A/L for forces F = .116, .197, .328piV, from bottom 
to top. x i s taken directly from experiment and r/ has been 
computed with A = 51.35nm and L — 15.61/im. The exper- 
imentalist's a is related to r\ by r\ — ujqAct = 27r-^- = 94.8a. 
The smaller points are the experimental results, the bigger 
points on the lowest curve are from Monte Carlo simulations, 
the full lines is the analytic study of the WLRC through the 
parametric representation (^|) with the values C/A given in 
the text. 

The regularization in the discrete model comes from 
the fact that the angle <j) — <t>' is defined modulo 2ir. 
The computation of < Xw > can be done explicitely 
at x — 0, and one finds that the continuum expres- 
sion L/A((l — cos0)/(l + cos#)} must be substituted by 
L/A((\ - cos0)/(l + cos0)i?(sin 2 9A/b)}, where the reg- 
ularization function is given in terms of Bessel functions 
by: R(x) = I\{x)/ Io(x). It is reasonable to assume that 
the discrete model is well approximated by the Hamil- 
tonian H r which is obtained from H by this same sub- 
stitution. We have computed the elongation properties 
of the WLRC from the ground state energy of the cor- 
responding regularized rod hamiltonian H r . The vari- 
ations of ^j^- with x now scale as a function of x/L, 
as in experiments. We introduce the intensive linking 
variable r\ — x A/L (related to the experimentalists' 
o by u — rj/ojoA). The partition function in (0) can 
be computed by the saddle point method in the limit 
L/i> 1 with 77 kept fixed. The saddle point is imagi- 
nary, k = in(a), and from its value one easily deduces the 
elongation of the chain using the general formulas given 
above. We obtain in this way the following parametric 
representation of the bell shape curves giving ^ versus 
77, for a fixed value of the force a. 
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An alternative to solving the wave equation associated 
with the regularized hamiltonian H r is to work directly 
with the discretized version of the effective energy (Q) , us- 
ing standard transfer matrix methods. We have checked 
that these two procedures give results for eo(a, — n 2 ) 
which are in very good agreement. Using interpolation 
techniques, one can eliminate k in M ) in order to ob- 
tain the ratio C/A in terms of -q, (z) / L and b/A. For 
various values of the force and of b/A, we have com- 
puted the empirical value of C/A obtained from each 
experimental point 77, (z)/L. With b/A = .12 all these 
values cluster nicely, allowing for a rather precise de- 
termination of C/A. For three values of the force, 
F = .116 pN; .197 pN; .328 pN, we find respectively 
C/A = 1.67±.12 ; 1.66 ±.10; 1.71 ± .09 (in this analysis 
we have restricted to the range of small enough super- 
coiling: \n\ < 1.5 for F = .197, .328 pN and \f]\ < 2.2 for 
F = .116 pN ). For each force, the value of C/A is the re- 
sult of a statistical weighted average and the quoted error 
is just the root mean square deviation, read off directly 
from the data. These three results agree and indicate a 
value of C/A ~ 1.68. However it should be stressed that 
the experimental points used here are preliminary and 
the study of systematic errors in the experiment is not 
yet completed so this value of C/A is not yet definitive. 
The bell shape curves, computed for each force with the 
corresponding value of C /A, are compared to the experi- 
mental points in figl. The agreement looks very satisfac- 
tory, indicating an overal consistency of the procedure. 
We have performed the same analysis with b/A = .05, .2. 
With b/A = .2 the agreement remains fairly good, but 
the value C/A for the three values of the forces exhibit 
fluctuations around 1.68 which can reach 20%. The value 
b/A = .05 appears to be excluded. 

Another possible method to get the ratio C/A could 
be to measure the curvature of the bell shape curve , 
r ?) = d 2 (< z > /L)/dr] 2 , at its maximum, which can be 
related to a\{a) — 2 (gjTeo(k 2 , a)) k Q through pertur- 
bation theory. The A/C is given in terms of the known 
function a\(a) M as: 
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At the moment the measurements of r j; is however 
not precise enough for a good determination with this 
method. 

We have shown that the WLRC must be regularised 
at small length scale. The corresponding model can be 
solved analytically and it accounts well for experimen- 
tal results at small supercoiling, giving a good method 
to determine the elastic constants ratio C/A. Obviously 
our theory is limited to the small force - small supercoil- 
ing regime. For instance the experiments show that for 



F > AbpN the extension is not symmetric for \ ~~ * — X- 
This kind of effect is totally beyond our simple elastic 
model which is intrinsically symmetric. Extending it re- 
quires the introduction of self avoidance to treat properly 
the plectoneme formation. The self avoidance will also 
naturally provide a cutoff length scale. However taking 
it into account properly is a major challenge 

We wish to thank J.-F. Allemand, D. Bcnsimon, V. 
Croquette and T.R. Strick for numerous exchanges, as 
well as J. -P. Bouchaud, A. Comtet and C. Monthus for 
useful discussions on the winding distribution of random 
walks. While completing this work, we learned that P. 
Nelson and his collaborators are working on the same 
problem in the limit of large forces. 
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